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Identification of a 12—17 day timescale in X-ray 
observations of GRS 1915+105 
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ABSTRACT 

Measurement of the energy emitted from accreting astrophysical systems provides an 
observational constraint on the plasma processes that may be operating within the 
disk. Here we examine the continual time variation over the past six years of the total 
X-ray flux from the microquasar GRS 1915+105. The application of differencing and 
rescaling techniques to RXTE/ASM data shows that the small amplitude fluctuations 
scale up to 12-17 days. A 17-day timescale in the X-ray fluctuations corresponds to 
half the measured binary orbital period of this system (33.5 ± 1.5 days). While this 
may be coincidental, it is possible that these two timescales may be linked by, for 
example, a turbulent cascade in the accretion disk driven by a tidally-induced two- 
armed spiral shock corotating with the binary system. Temporal scaling is found only 
in the ever-present small fluctuations, and not in the intermittent larger-amplitude 
fluctuations. This is consistent with the basic model for this source which consists of 
a steady, cold outer disk and an unstable inner disk. 
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1 INTRODUCTION 

GRS 1915+105 was discovered in 1992 by the WATCH All- 



the proportionality of dynamic timescales to mass implies 
that thousand-year events in a quasar could be studied over 



sky X-ray Monitor on board the GRAN AT satellite ( Castro- 



Tirado et al. 


1992 


), and lies in the Galactic plane at a dis- 


tancc of 11-12 kpc ( 


Mirabel and Rodriguez 


1994 


; Fender et 



minut es in a microquasar (see for example Mirabel et al 

(EH»- 



al. 



1999 ). The extreme variabili ty of the X-ray flux f r om it s 



199E 



accretion disk, first noticed by Castro-Tirado et al 
coup l ed- with its jet o mission (Mirabel and Rodriguo 

la 



hagac et. aj.20flCj-Chaty pt al l 2flfU), makes it a 



(1994) 



Spectral states, quasi-periodic oscillations (QPOs), and 
the rich temporal variability of GR S 1915+105 arc the focu s 
of much resea rch (s e e for example |Nayakshin ct al. ( 200C ) 



1094 



Belloni et al 
al 



fl200l|); [Klein- Wolt et al] Q2002j ); |Rodrigucz ct 



(2002); Ueda et al. (2002); and references therein). Be- 



ideal source f or the study of accr etion physics. At 14 + 4 

hereafter referred to as 



solar masses (Greiner et al 



2001 



G2001), it can thus be considered as a Galactic stellar-mass 
analog ue of a maooivo oxtragalactic black hole oyotom ouch 



fore the measurements by G2001, mass estimates for GRS 
1915+105 made use of the tem poral and spectral similarity 
to microquasar GRO J1655-40 ( |Grove et al. 1998), another 
Galactic X -ray source with superluminal radio jets (Zhang 



as a quasar, hence the term microquasar — The importance 

of microquas ars as accessible astrophysical laboratories is 
discussed by Mirabel and Rodriguez ( ]199E| ); for example, 



et al 



1994 ) together with a binary mass function appropri 
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ate for a black hole (Bailyn et al. 1995). The radial veloc- 
ity measurements performed by G2001 determined a mass 
of 14 ± 4 solar masses, leading the authors to suggest that 
the X-ray variability can be accounted for by instabilities 
in a radiation-pressure-dominated accretion disk near the 
Eddington limit. The observations of G2001 also constrain 
the mass of the donor star to 1.2 ± 0.2 solar masses, and its 
orbital period around the black hole to 33.5 ±1.5 days. 

There is now much evidence for the basic model of a 
cold, steady outer disk and an unstable, radiation-pressure- 
dominated inner disk in GRS 1915+105 (see for example 
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Belloni et al.| ([L997a|,|b[) and references therein). Hence there 



are two emitting regions to consider: an inner region that 
empties and refills on timescales of seconds , perhaps due 



to a rapid local viscous-thermal in stability (Lightman and 



Eardlejj |l974 



Belloni et al 



1997a Janiuk et al, 2000); and 



a steadier outer region varying only over the much longer 
timescale t v i ac associated with viscous dissipation of bulk 
plasma flows. We may therefore expect t V i SC to emerge as the 
longest timescale over which correlation may be present in 
observations of t he total lum inosity; a recent study suggests 
tvisc ~ 10 years jVilhvj |2002j) . 

In this paper we apply the technique of differencing and 
rescaling to the GRS 1915+105 X-ray data, and show that 
this implies the existence of a fundamental timescale for the 
system in the range 12-17 d ays. The te chni que of differenc- 
i ng an d rescaling was used by (She- (1991) and She and Orszag 
(1991) to stud y a physical model of in t ermit tency in turbu- 
lence, and by Mantegna and Stanley (1995) to investigate 
fluctuations in the value of a fi nancial index. Mor e recen tly, 
the method has been used by Greenhough et al. (2002) to 
compare ac c reting with non-accreting X-ray sources, and by 



Hnat et al 



(2002) to quantify turbulence in the solar wind. 



In summary, we first construct a set of differenced series 
Z (t, t) from the original series y(t), for a range of values of 
the time-lag r: 



Z(t,r)=y(t)-y(t-r). 



(1) 



A set of PDFs is then calculated for the amplitude of the 
differenced series Z, one for each value of r, denoted by 
P(Z,t). The best-fit slope of logP(0,r) against logr de- 
fines the scaling exponent m. We then seek a common func- 
tional form of these PDFs by rescaling both axes such that 



Zr m = Z s 



and 



Pr~ m = P s 



(2) 



If the separate PDFs for values of r up to some T ma x collapse 
onto one curve, this implies that the X-ray fluctuations are 
controlled by a single physical process on timescales below 

This technique enables us to test in particular for 
t he joint presence of temporal scaling and non-Gaussianity 



( Hastings and Sugihara 



Sornette 



1993 



Kantz and Schreiber 



1997 



2000), which together are st rong indications of cor- 



related processes such as turbulence (Bohr et al. 199S ) . GRS 
1915+105 is ideal for this type of statistical analysis in that 
(i) it remains active on very long time scales, and (ii) this 
time scale is very well separated from those associated with 
the inner disk. 



2 SCALING OF GRS 1915+105 X-RAY DATA 
2.1 Summary of the data 

The raw X-ray data, provided by t he All-Sky Monito r 
(ASM) on board the RXTE satellite ( jSwank et al~| |200l| ), 
are held at the Goddard Space Flight Center (GSFC) and 
can be accessed via their website^. Calibration is undertaken 
by the ASM/RXTE team a nd the processed d ata are freely 
accessible on their websiten (|Bradt et al. 2001), where there 



1 http:/ /heasarc. gsfc.nasa.gov/docs/xte/asm_products. html 

2 http://xte.mit.edu/XTE/asmlc/ASM.html 
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Figure 1. X-ray time-series for GRS 1915+105, 1996 Feb. 20 
- 2002 Apr. 12 (MJD 50133-52376); mean sampling interval 94 
minutes. There are 32,800 data points, each of which is the num- 
ber of counts in the range 1.3—12.2 keV obse rved in a 90-second 



interval. Data ob tained by ASM on RXTE (Swank et al. 
Bradt ct al [2001 ). 



2001 



are 32,800 observations over the period 1996 February 20 
to 2002 April 12, which we plot in Fig. [j]. Each point rep- 
resents the total X-ray flux (measured by the number of 
counts during periods that last 90 seconds) integrated over 
the range 1.3-12.2 keV. The breakdown into three energy 
bands (1.3-3, 3-5 and 5-12.2 keV) is also available; we have 
carried out separate analyses of these components, which 
yield very similar results to those for the total flux which we 
present here. The source luminosity is sufficiently high to 
permit the neglect of instrument thresholds, uncertainties, 
and other sources in the field of view. Sampling intervals be- 
tween the 90-second X-ray counting periods are distributed 
with a mean of 94 minutes, and 90% of the intervals are 
below 193 minutes; the exact observation times are known 
and therefore we do not assume even sampling in our anal- 
ysis. When differencing the X-ray intensities to construct 
P(Z,t), we exploit the fact that each measured intensity is 
obtained at a precise time that is stated in the dataset to 
within 0.1s. We use these sampling times to find intensities 
that were observed within 2 ±01 r of the chosen value of r. 
While the ASM data is not filtered prior to sampling, and 
hence the true frequency spectrum is distorted by aliasing, 
the differencing technique uses only time domain informa- 
tion that accurately represents the true signal at the given 
times. We do not attempt any kind of interpolation. 



2.2 Results 

In Fig. ^(a) we plot the probability of the smallest fluctua- 
tions, P(0,r), as a function of r. The straight lines on the 
double logarithmic axes indicate that P(0, r) is a power- 
law function of r up to r ~ 2 7,5 . This corresponds to a 
timescale of 12 days, since this value of r represents 2 7,5 
times the mean sampling interval of 94 minutes. This break 
in scaling appears to be a robust feature of the data. For 
example, Fig. ^|(b) shows the variance a 2 (r) of the set of all 
fluctuations for a given value of r; that is, the variance of 
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Figure 2. (a) Dependence of the probability of the smallest- 
amplitude fluctuations (differences), P(0, r), on differencing time 
interval r. This is a logarithmic plot of P(0, r) against r = j'tq, 
where j steps up in integer powers of two and to = 94 minutes. 
Vertical error bars indicate a factor of ± 1/ \/~N, where N = num- 
ber of pairs of data points separated by r that are used to cal- 
culate P(0,t); horizontal error bars of ±2 01 indicate the range 
of t separating data points used to obtain P(Z,t). The slope 
m = —0.08 is used to derive Fig. 3. (b) Dependence of the vari- 
ance, ct 2 (t), of the set of differences (Z,r) on differencing time 
interval r; error bars are calculated as in (a). The slopes of the 
solid lines in both plots are calculated from linear regression with 
95% confidence intervals. Taken together, the plots suggest that 
scaling extends from 0.5 days to f2-17 days. 



the differenced data from which the PDFs P(Z, r) are cal- 
culated. It implies the existence of scaling for r up to 2 s , 
corresponding to 17 days. Our key result is that there is a 
clear break in scaling between the 7.5th and 8th octaves, 
corresponding to a timescale of 12-17 days. We emphasise 
that 12-17 days is not a periodicity in the X-ray luminos- 
ity but a timescale at which a break in the scaling of its 
statistical properties is observed. In Figs. §(a,b) a second 
scaling region emerges between octaves 7.5 and 10.5. This 
represents a rather restricted range of timescales, however, 
and will not be considered further in this paper. 

We next use the slope m = —0.08 of the linear fit over 
the lower range of timescales in Fig. ^(a) to carry out rescal- 




Figure 3. (a) Rescaled PDFs of differenced time-series, where 
Z a denotes the differences rescaled according to r and m as at 
Eq. |^. Curves are plotted for time intervals r = jro where j steps 
up in integer powers of two and to = 94 minutes; j m ax = 8 (= 17 
days). The PDFs collapse onto one non-Gaussian curve between 
the dot-dash lines at | Z B [= 30, indicating that only the small 
fluctuations exhibit temporal scaling. For comparison purposes, 
the dashed line is the Gaussian fit to curve j = 1 for all Z s . 
(b) Peak region of (a) enlarged to display detail, with Gaussian 
(dashed line) fitted to curve j = 1 for | Z s \ < 30 only. 



ing, as described at Eq. g, for fluctuations of size Z. Fig- 
ure ^a) shows a collapse onto one curve for | Z s |< 30, im- 
plying that small fluctuations follow the same scaling, but 
large fluctuations beyond | Z s |~ 30 clearly do not scale 
in this way. In Fig. @(b) we show that the scaling regime 
indicated in Fig. pi a) is peaked, long-tailed, and definitely 
non-Gaussian, which is suggestive of turbulent behaviour 



(Bohr et al 



1995 



We now show that the small fluctuations that exhibit 
scaling are continually present, whereas the large fluctua- 
tions are not. Figure provides representations of: (a) the 
full X-ray time-series for GRS 1915+105; (b) the small fluc- 
tuations that lie between the dot-dash lines (| Z a |< 30) in 
Fig. [| and (c) the large fluctuations (| Z s |> 30) only. Dots 
are the data points at which the fluctuations terminate, that 
is the latter points of the pairs that are differenced to con- 
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Figure 4. X-ray tinre=series plullud using duly Lu identify the 

origin of fiuctuations of different sizes; (a) aii fiuctuations (all 
data points), (b) small fluctuations only (between the dot-dash 
lines in Fig. 3), (c) large fluctuations only (outside the dot-dash 
lines in Fig. 3). This plot shows that only the small fluctuations 
originate from all parts of the time-series, and Fig. ^ shows that 
only these fluctuations exhibit scaling. 



struct P(Z,t). Thus, for example, in Fig. ||(a) there is a 
dot for every raw data point, whereas in (b) we plot one 
dot only for each pair of intensities whose rescaled differ- 
ence Z s lies between Z s = ±30. It is clear that the small 
scaling fluctuations arise from all parts of the time-series, 
whereas the larger non-scaling fluctuations are sometimes 
absent (notably for long periods in 1997 and 1998-99). 



3 DISCUSSION AND CONCLUSIONS 

Using the differencing and rescaling technique, 12-17 days 
emerges as a resilient timescale in the total X-ray luminosity 
(1.5-12 keV) of GRS 1915+105, and is shown to originate 
in only the small luminosity fluctuations. This is consistent 
with the accepted model of a cold (low emission and vari- 
ability) outer disk that is stable over (viscous) timescales 
well above those of the inner disk; only in the outer disk 
could processes evolve without interruption over periods of 
days. 

It is also intriguing that 17 days is approximately half 
the binary period (33.5 ± 1.5 days) of this source. Although 
this may be coincidental, it is interesting to speculate on its 
possible physical significance. For example, this result lends 
some support to an interpretati on in terms of tidally -driven 
waves which we now describe. 



Lubow et al 



(1993) 



fluid experiments beginning in the 1970s, in which it was 
found that a dynamical ins tability can arise i n a fluid whose 
base state is non-circular; |Goodma"n (1993) originally ap- 
plied this idea to tidally-distorted accretion disks in close 
binary systems. The instability is caused by a perturbation 
in the frame rotating with the binary at angular frequency 
Q b ■ Trailing waves can couple via the tidal mode to generate 
a leading wave, which can then reinforce the original distur- 
bance via the tidal distortion. Vishniac and Zhang (199£) 



ings, and Ryu et al. (1996) show how the instability can 



bec ome global (for oth e r typ es of instability see, for exam- 
ple, |Caunt and Taggerj ( J200l| )). 

In the mode with the largest amplitude, the surface den- 
sity varies slowly with radius and has an angle and time de- 
pendence given by sin(2#— 2f2si) (Lubow et al. 1993). Differ- 



ential rotation within the disk shears this disturbance into 
a two-armed spiral shock that co-rotates with the binary; 
amplification caused by the inwardly -increasing velocity is 
balanced by dissipation ( Blondin 200C). This tidally- induced 
shock wave has been shown numerically to b e a robust fea- 



ture of the accretion disks of binary systems ( Blondin 



200C 



and references therein), and their possible contribution to 
the light curves of in termediate polars was demonstrated by 
Murray et al. (1999). There exist, in addition, ma ny obser- 



vation s of sp iral shocks in binaries (see for example Baptista 



et al 



-f 

(2003) and references therein) and in AGN (Rouan et 



al.|(|l998D~ and references therein, for example). It is plausible 
that the periodic stirring effect of the rotating shock could 



drive a turbulent cascade ( Soodman 



1993 



Bohr et al 



199£ 



in which plasma is heated by viscous dissipation within vor- 
tices covering a range of spatial and temporal scales smaller 
than those of the driving process, and with non-Gaussian 
fluctuations. It is this dissipation, spatially integrated across 
the accretion disk, that drives the X-ray emission. Our ob- 
servations are consistent with the hypothesis that in the 
present case of turbulence in a bounded system, the tempo- 
ral scaling properties - and in particular the upper timescale 
of 12-17 days - persist in the X-ray signal. With two spiral 
arms each rotating with a period of Tb = 2n/Qs, one would 
therefore expect to observe temporal scaling up to Tb/2 in 
observations of the total luminosity. For GRS 1915+105, 
T B /2 = (33.5 ± 1.5) /2 = 17+1 days. 

Our principal result thus provides quantitative con- 
straints on models for transport and heating in the accretion 
disk of GRS 1915+105, and our analysis also suggests pos- 
sible physical mechanisms that may be at work. 
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